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Abstract. We consider 2p ^ 4 order differential operator on the real line with a periodic 
coefficients. The spectrum of this operator is absolutely continuous and is a union of spectral 
bands separated by gaps. We define the Lyapunov function, which is analytic on a p-sheeted 
Riemann surface. The Lyapunov function has real or complex branch points. We prove the 
following results: (1) The spectrum at high energy has multiplicity two. (2) Endpoints of all 
gaps are periodic (or anti-periodic) eigenvalues or real branch points. (3) The spectrum of 
operator has an infinite number of open gaps and there exists only a finite number of non-real 
branch points for some specific coefficients (the generic case). (4) The asymptotics of the 
periodic, anti-periodic spectrum and branch points are determined at high energy. 

1. Introduction and main results 
Consider the self-adjoint periodic operator H acting in L 2 (R) and given by 

d 2p ^ d j d j 

H = H + q, H = (-iy—, g = £ P> 2, (1.1) 

3=0 

q 3 eL l real {T), j e% = {l,2,...,p}, T = R/Z, (1.2) 
where Z is the set of all integers. Let Wj(R),j G N = Z fl [l,oo), be the Sobolev space of 

functions /, G L 2 (R). Here and below we use the notation f — %, = We define 
the self-adjoint operator H using the quadratic form with the form domain Dom fd{H) = 



W 2 (R) (see Proposition IQh 

It is well known (see |DSj . Ch. XIII. 7. 64) that the spectrum cr{H) of H for the sufficiently 
smooth coefficients qj, j G N p , is absolutely continuous and consists of non-degenerated inter- 
vals & n ,n = 1, ...,N G ^ oo. These intervals & n and & n+ i are separated by the gap g n with 
length |fj n | > and Nq — 1 is a number of the gaps. Theorem 11.11 extends this result to the 
larger case qj G L 1 (T). 

The typical applications of our operator H are the vibrations of beams, plates and shells: 

(1) The standard Kirchhoff-Love model of the bend of beams and plates provides the Euler- 
Bernoulli equation y"" = Xay (see [TY Wj . Ch. 5.9). 

(2) The Vlasov model of the bend of cylinder shells (see [NCMj . Ch. 1.1.14) gives the 
equations of vibration having the form + b\y = Xby. 

Here y is the normal displacement of the plate (or shell), the functions a (or b, b\) are defined 
by the parameters of the plate (or shell): Young's modulus, Poisson's modulus, rigidity and 
thickness. 
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The high order differential operators arise in the inverse problem method of integration of 
non-linear evolution equations. There exist the Lax pairs, where the self-adjoint operator is a 
high order operator and the corresponding non-linear Lax equation is integrable by the inverse 
problem method, see jDKNj . Many physically interesting equations have this form, see |AC] . 

Recall that the spectral theory for the high order operators with decreasing coefficients is 
well developed, see [BDT], [Su] and the references therein. The results for high order periodic 
operators are still modest. 

We describe our goal. In the case p — 1 the spectrum of the Hill operator — jm + Qi is a 
union of spectral bands, where all endpoints of the bands are 2-periodic eigenvalues of the 
equation — y" + q\y = \y. In the case p = 2 the spectrum of the operator H is also a union 
of spectral bands, but all endpoints of the bands are 2-periodic eigenvalues of the equation 
yim + qy = \y or the branch points of the Lyapunov function [BK2J. Until now there are 
no any results about the multiplicity of the spectrum at high energy, number of gaps (is it 
finite or infinite ?), asymptotics and type of endpoints of the gaps at high energy etc for the 
operators H, p > 2. Our main goal is to answer some of these questions. 

In order to describe our results we consider the equation 



-l)PyW +qy = \y, q 



j=0 



(f,A)e 



x 



;i.s) 



where C is the complex plane. If all coefficients q^q^- ^ € L 1 (T), then the standard mon- 
odromy matrix is well defined (see |DSj . Ch. XIII. 7). If some coefficient qj G -^(T), q'j ^ ^ 1 (T), 
then the standard monodromy matrix is not well-defined, since, in general, the derivative of 
y(2p-i) j s no t continuous. In this case we will introduce the modified symplectic monodromy 
matrix, see (jl.iup . We think that it will be convenient even for smooth coefficients qj. We 
rewrite the equation (II. 3p in the vector form by 



Y' - V(X)Y = QY, (t, A) € K x C, 
sec [Naj. Ch. II. 4, where the vector- valued function Y is given by 



;i.4) 
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and the 2p x 2p matrices- valued functions V, Q are given by 
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O m) „ is the m x n zero matrix, l n is the n x n identity matrix. If all qj G L 1 (T), then 
Q G L 1 (T) and there exists a 2p x 2p matrix- valued solution Ai(t, A) of equation (11. 4p with 
the initial condition A4(0, A) = l 2p . In this case the modified monodromy matrix A4(l, A) is 
well-defined and entire. Its characteristic polynomial D is given by 

D(t, A) = det(M(l, A) - rl 2p ), (r, A) G C 2 . (1.7) 

An eigenvalue of -M(l, A) is called a multiplier, it is a zero of the algebraic equation D(-, A) = 0. 
Each A4(1,A),A G C, has exactly 2p (counted with multiplicities) multipliers Tj(X),j G N 2p . 
Due to (11.101) . the matrix M. is symplectic. Then r is a multiplier iff r _1 is a multiplier. The 
multipliers have asymptotics 

Tj (X) = e z ^(l + 0(\z\~ 1 )) as |A| ->• oo, A G C+ = {A : ImA > 0}, all j G N 2p , (1.8) 
argA G 



A^ G 5 



zGC: are; z G ( — 



2p 2p 



see Lemma [4.11 Here and below Uj are the zeros of the polynomial u 2p — (— l) p labeling by 

W 2 j-1 = — L02(p-j+l) , = — ^2p-j+l, j € N p . (1.9) 



—i — i 

e p . 



el- 



even p 
e ll 2p^~ l \ odd p 



Note that w p+ i = — u p = i and if p is even, then uj\ = —uji v = 1. 

The coefficients of the polynomial D(-,X) are entire functions in A. It is well known (see, 
e.g., [FoJ, Ch. 8) that the roots Tj(X),j G N 2p , constitute one or several branches of N F ^ 1 
analytic functions that have only algebraic singularities in C. Asymptotics ( II. 8p show that 
Np = 1, i.e. Tj are branches of the unique function r analytic on the 2p sheeted Riemann 
surface. Moreover, these asymptotics define the functions Tj in C + for |A| large enough. The 
detailed results about the branches Tj will be given in Section 4. 

We formulate our first preliminary results. 

Theorem 1.1. The monodromy matrix M. = -M(l, •) is symplectic, i.e. it satisfies the identity 



M T JM = J, where J 



1)% 



J, 
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[1.10) 



and J~ T = — J . Furthermore, there exists an analytic junction A on the connected p-sheeted 
Riemann surface & having the following properties: 
i) All branches of A have the form Aj 

$(v,X) = f[(v-A j (X)), v= T -± T — } (r,A)GC 2 , r^O, (1.11) 
j'=i 



(2r)f 



\{jj + Ta G N p , and satisfy: 



Aj(X) = coshzcoj + O 











as 









X — > oo, A G 



;i.i2) 



ii) If Aj(X) G (—1,1) for some (j, A) G N p x M, and X is not a branch point of Aj, then 

a;.(a)^o. 

Hi) The spectrum cr(H) of the operator H satisfies 

a(H) = a ac (H) = {X G R : Aj(X) G [-1, 1] for some j G N p }. (1.13) 
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FIGURE 1 . The Riemann surface of the Lyapunov function for p = 3 for the case of small 
coefficients. Identities (l4~22|) give A 2 (r^„) = A 3 {rf n ) and A^rf j = A 2 (r^J. Then the 
second and third sheets of the surface are attached along the cuts (rf n _ x , rjf n ), n 6 N. We 
attach the upper (lower) edge of each cut on the second sheet to the lower (upper) edge of 
the same cut on the third sheet. Similarly, first and second sheets of the surface are attached 
along the cuts (r^" n _ 1 , rj n ), n G N. Thus, whenever we cross the cut, we pass from one sheet 
to another. 

Remark. 1) If A G o'(H), then some branch Aj(A) is real and the corresponding multiplier 
Tj(X) is complex and |tj(A)| = 1. It is more convenient to study the real function Aj(A) on 
the spectrum <j(H), than the complex multiplier r 3 -(A) on cr(H). 

2) The surface Si is connected. For the first and second order operators with the matrix- 
valued potentials the corresponding surface may be disconnected (see |CKj . |K1] . |K2] ) . 

3) The proof of i), ii) repeats essentially the argument from [CKj . |K1] . [K2J. 

4) The monodromy matrix for the second order operators has asymptotics in terms of 
cos and sin bounded on the real line. The monodromy matrix for high order operators has 
asymptotics in terms of cosh and sinh, see (13.111) . unbounded on the real line. 

The zeros of D(l,-) (or D(— 1, •)) are periodic (or antiperiodic) eigenvalues of the equa- 
tion (— l) p y( 2p ^ + qy = Xy, where y are 1-periodic (or 1-antiperiodic) functions. Denote by 
XQ,Xf n ,n ^ 1, the periodic eigenvalues and by Xf n _ 1: n ^ 1, the antiperiodic eigenvalues 
labeling by (counted with multiplicity) 

Aj ^ A 2 < X$ ^ A 4 < Xt ^ A 6 ^ X± ^ Xf ^ A3 ^ X$ ^ A5 ^ A^ .... 

For the polynomial $ given by (11.111) we introduce the discriminant p(X), X G C, by 

p= J] (A^.-A,) 2 . (1.14) 

A zero of p is a ramification point (or simply a ramification) of the Lyapunov function A. 

Remark. 1) Ramification is a geometric term used for 'branching out', in the way that 
the square root function, for complex numbers, can be seen to have two branches differing 
in sign. We also use it from the opposite perspective (branches coming together) as when a 
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covering map degenerates at a point of a space, with some collapsing together of the fibers of 
the mapping. 

2) Recall that all endpoints of gaps of the spectrum of the Hill operator (i.e., p = 1) are 
periodic or anti-periodic eigenvalues. The situation is more complicated for the high order 
periodic operators and the periodic operators with the matrix potentials. In these cases 
the endpoints of gaps are periodic or anti-periodic eigenvalues, or ramifications (zeros of the 
function p). The numerical analysis for the fourth order operators and the second order 
periodic operators with the 2x2 matrix potential shows that ramifications can be non-real 
for some values of the coefficients and can become real and to create the gap for some other 
values of the coefficients. This behavior is similar to the behavior of the resonances in the 
scattering problem for the Schrodinger operator (see, e.g., [K5j . [Z]). In fact, this was a reason 
for us to use the term resonance for the zero of the function p for the periodic operators in 
our previous papers [BK1] . |BK2j . [BBK]. |CK] . ... But now we will use the term ramification 
for such points because the term resonance is overloaded and is used in different other senses. 

We shortly describe the unperturbed operator H° = (— see more in Section 2. The 
unperturbed multipliers t®, the Lyapunov function A with all branches A° are given by 

r °(A) = e^ z , A°(A) = cos A^, A°(A) = cosh zu v (j, A) G N 2p x C+. 
The unperturbed discriminant p° has the form 



P° 



- YJ_ (cosh zuj — cosh zujh) 2 . (1.15) 



The 2-periodic eigenvalues A^ = (irn) 2p ,n ^ 1, have multiplicity 2 and the periodic eigen- 
value Aq' + = has multiplicity 1. The function p° is entire and has the zeros r® n , k G N p _i, n 
0, given by 

rln= - 1 — > c fc = cos— , 1 > ci > c 2 > ... > c p _i > 0. 1.16 
\c k / 2p 

The zero A = of the function p° has the multiplicity p — 1 and each another zero has the 

multiplicity 2. The spectrum o~(H°) has the multiplicity 2. The Riemann surface &° of the 

Lyapunov function A for the operator H° coincides with the Riemann surface of the function 
i 

\p with the unique branch point at A = 0. 
We determine the sharp asymptotics of the ramifications. 

Theorem 1.2. i) The function p is entire, real on R and satisfies 

p(\) = p°(\)(l + 0(\z\- 1 )) as|A|^oo, \X-r° kn \>l, V (k, n) G x N. (1.17) 



ii) The function p has the zeros r£ ,r kn , (k,n) G N p _i x N, which satisfy 



rt n = (-l) k [-\ 



(l + ^ [(-lT + %,0 ± c k \q P ,n\ + O(i)]) (1.18) 



as n — > oo, where 

»i 

-i2nnt 



q P ,n= / q P (t)e- Umt dt, n^O. (1.19) 
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Remark. 1) Asymptotics (I1.17P show that p ^ 0, since p° 7^ 0. Note that for the second 
(and first) order operators with the p x p matrix-valued potential the corresponding function 
may be equal to (see |CK], [Kl], |K2] ). 

2) The Riemann surface roughly speaking, is close to at high energy. In general, the 
points r^ n are simple branch points of A (square root type) for n large enough. The surface 
for p = 3 is shown by Fig. [TJ 

We describe the structure of the bands and the gaps at high energy. 

Theorem 1.3. i) The branch A p is real analytic function on the interval (Ai,oo) for some 
n G N and A+ < A~ o+1 ^ A+ o+1 < K a +2 < K +2 < ■■■■ Moreover, if n ^ n , then each 
interval [A+, A~ +1 ] is a spectral band with the spectrum of the multiplicity 2, and each interval 
(A~, A+) 7^ is a gap and A 2 (A^) = 1. There are no other bands of H to the right of A+ Q . 
ii) The periodic and antiperiodic eigenvalues A^ satisfy: 

Xn = (vrn) 2p (\ + (-l) p+1 g P ,o ± \%,n\ + ^ ) as n^oo. (1.20) 

Remark. 1) The spectrum of H has multiplicity 2 at high energy. The spectrum of the 
Schrodinger operator with the p x p matrix- valued potential and the first order operator with 
the 2p x 2p matrix- valued potential has multiplicity 2p at high energy (see |CK] , |K1] ). 

2) The periodic and antiperiodic eigenvalues accumulate at +00. The ramifications accumu- 
late at ±00. But for second (or first) order systems the periodic and antiperiodic eigenvalues 
and the ramifications accumulate at +00 (or ±00), see |CKj . |K2j (or |Klj ). 

3) The spectrum of H at high energy is described by the branch A p of the Lyapunov function. 
The structure of the spectrum at high energy is similar to the structure of the spectrum of 
the Hill operator: 

(a) the spectra of H and of the Hill operator are similar as the sets, including multiplicities; 

(b) endpoints of gaps are periodic and antiperiodic eigenvalues only; 

(c) the sharp asymptotics of these eigenvalues are expressed in terms of the Fourier coeffi- 
cients of the potential. 

Recall that there exists an infinite number of open gaps in the spectrum of the first and 
the second order operators for some specific potentials |CK] , |Klj , |K2j . Now we describe this 
situation for our case. 

Corollary 1.4. Let the coefficient q p satisfy \ qp(t)e _,27mjt <it| ^ for some infinite se- 
quence of indices nj — > 00 and some < a < 1 . Then 

i) All gaps j nj = (A^, A+) 7^ are open and the gap-length \j n:j \ — > 00 as j — >■ 00. 

ii) All ramifications rjjr n ., k G N p _i are real and there exist an infinite number of the non- 
empty intervals (r^ n ., rj£ n .) C K with the length \r^ n . — r^ n .| — > 00 as j — > 00. 

Note that if | q v {€)e~ l2l * nt dt\ ^ as n — > 00 and for some < a < 1 (generic periodic 
coefficients). Then there exists only a finite number of non-real ramifications (branch points 
of the Lyapunov function) and all high energy gaps are open. 

A great number of papers is devoted to the inverse spectral theory for the Hill operator: 
Dubrovin jD], Garnett and Trubowitz |GT] , Its and Matveev |IM] , Kappeler |Kap| , Kargaev 
and Korotyaev |KK] . Korotyaev |K3j . Marchenko and Ostrovski |MO] , Novikov |No] etc. Note 
that Korotyaev [K4] extended the results of |MO] , |GTj . [Kap],[KK], |K3] for the case —y" + qy 
to the case of periodic distributions, i.e. — y" + q'y on L 2 (IR), where periodic q G Lf oc ( 
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We describe now the results for vector differential equations. Recently the inverse prob- 
lem for vector-valued Sturm-Liouville operators on the unit interval with Dirichlet boundary 
conditions, including characterization, was solved by Chelkak, Korotyaev |CK1] . |CK2j . The 
periodic case is more complicated and a lot of papers are devoted only to the direct problem 
of periodic systems: Carlson |Calj . |Ca2j . Gelfand and Lidskii |GL] . Gesztesy and coauthors 
|GLj . Korotyaev and coauthors [CKJ, [BBK], |Klj . |K2j . etc. We describe results for first and 
second order operators with the periodic p x p matrix- valued potential from [UK], |Klj . |K2j : 

1) the properties of the Lyapunov function, defined on the Riemann surface, are described, 

2) the conformal mapping with real part given by the integrated density of states and 
imaginary part given by the Lyapunov exponent is constructed and the main properties are 
obtained, 

3) trace formulas (similar to the case of the Hill operators) are determined, 

4) an estimate of gap lengths in terms of potentials is obtained, 

5) sharp asymptotics of periodic eigenvalues and ramifications are determined. 

Note that the discrete periodic systems were studied in |KKulj . |KKu2j . The results for first 
and second order operators are important for us, since we plan to repeat one for even order 
periodic operators. In fact this is the motivation of our paper. Note that the case of even 
order periodic operators is more complicated than the case of first and second order operators, 
since in the first case only one fundamental solution is bounded on the real line and all other 
fundamental solutions are unbounded on the real line. 

We describe the fourth order operators H = <9 4 + dq 2 d + <?i- The results for decreasing 
coefficients are more developed, see |AP] . |GM] . |HLO] . |LO] . We mention the paper [CPS], 
|McL] about the inverse problem for fourth order operators on the unit interval. Now we 
describe the periodic case. The authors |BK2] obtained the following results for the operator 
H = <9 4 + dq 2 d + qi (the case q 2 = see in |BK1] ): 

(1) The properties of the Lyapunov function on the 2-sheeted Riemann surface are described. 
The asymptotics of the spectral gaps and ramifications are determined at high energy. 

(2) If q 2 = 0, qi — > or q\ = 0, q 2 -» 0, then there exists a small non-empty spectral band 
with the spectrum of multiplicity 4. The beginner of this band is the ramification, which 
coincides with the top of the spectrum. The spectrum in all other bands has multiplicity 2. 

(3) There exist both real and non-real ramifications for some specific potentials. 

The spectral properties of the periodic Euler- Bernoulli equation (ay")" = Xby were studied 
by Papanicolaou |Plj . |P2j . |PK] (jointly with Kravarritis). It was shown that the spectrum 
is a union of non-overlapping bands of multiplicity 2, similar to the case of the scalar Hill 
operator. The beginning of the spectrum is both a simple periodic eigenvalue and a branch 
point of the Lyapunov function. All other ramifications are negative. 

Consider the operator H with p ^ 2. The old well known results see in the book [NaJ. 
Tkachenko |Tk] obtained the eigenfunction expansion formula for the operator H. Mikhailets 
and Molyboga [MMlJ, |MM2] determined asymptotics of eigenvalues for the operator (— l) p d 2p + 
q on the circle T = M./Z, where q is a distribution. Galunov and Oleinik |GO] considered the 
operator (— l) p d 2p + 5 per on the real line, where 5 per is a periodic 5-function. 

It is important that for p = 1 the spectral analysis of the operator on the circle (the periodic 
and antiperiodic spectrum) is equivalent to one of the operator H on the real line. The main 
tool is the analysis of the entire Lyapunov function. The situation for p ^ 2 is much more 
complicated (see [BKlj . |BK2j ). In this case the Lyapunov function A has the complicated p 
sheeted Riemann surface. 
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FIGURE 2. The numbers u)k for p — 2 and p — 3. The domains {zuJk, ^ argz ^ ^} are shaded. 

In the present paper we extend some of results from [BKlj . [BK 2j about the case p = 2 
to the case p ^ 2. We construct the Riemann surface for the Lyapunov function of H and 
describe this surface for large |A|. Moreover, we determine asymptotics of the ramifications 
and periodic and antiperiodic eigenvalues at high energy. 

The plan of the paper is as follows. In Sect. 2 we describe the multipliers for the unperturbed 
operator. In Sect. 3 we describe the basic properties of the monodromy matrix Ai. In order 
to determine the asymptotics of the monodromy matrix Ai at high energy we use so-called 
Jost type solutions with "good" asymptotics at high energy. In Sect. 4 we obtain the main 
properties of the multipliers, the Lyapunov function and the function p and prove Theorem 
11.11 Moreover, we consider some simple examples. In Sect. 5 we prove our main Theorems II. 21 
and 11.31 In the proof using the mix of arguments both for the fourth order operator |BKlj . 
[BK2J and for the systems |CKj . |Klj . |K2j . we determine the asymptotics of the ramifications 
and periodic eigenvalues analyzing directly the determinant D of the monodromy matrix in 
the neighborhoods of ramifications (see Lemma T5.2I and the proof of Theorem ll.2[ 11.3]) . In the 
end of Sect. 5 we prove the simple Corollary 11.41 from Theorem 11.31 Some technical proofs 
are placed in Appendix. 

2. Properties of the unperturbed operator 
The numbers Uj, given by (11. 9p . satisfy 

p odd : Recoup = Reco>2p-i < ••• < Reo;4 = Reu^ < Reuj2 = Reui, 

lmw 2 j-i<0, uj 2 j = w 2 j-i, all j G N p , (2.1) 

p even : —1 = U2 P < Reu2 P -i = Reo;2p-2 < ••• < Rec^s = Reo;4 < Reo;3 = Reo;2 < toi = 1, 

lmui2j < 0, UJ2 j+ i = uJ2j, all j G N p _i, (2.2) 

see FigfSJ Moreover, 

Up+j+i = erfj, u p+j = EjTij, u p+j - uj p+j+l = (-ly+^iCjfjj, (2.3) 
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for all j = —p + 1, — p + 2, p — 1, where 

7T7 f ie^p , j even j 1, j even 

2p y-ie 2 p, j odd [e 2p , j odd 

We introduce the step functions f^(A), (j, A) G N 2p x C, which are constant in each half-plane 
C± = {A G C : ± Im A > 0} and given by 

n i (A) = |_ Wi ' ImA ^° , (2-5) 

I uij = uj£ } Im A < 

where 

f ' + 1 ' odd \ 3 ' + 1} J eVen ' J ^ 2P 

odd p : £ = < ' , even p : I — < j — 1, j odd, j 7^ 1 . (2.6) 

I J - 1, J even . . 

k ( J, J = Uj' 

Lemma 2.1. T/ie functions Oj(A), A G C, j G N 2p , satisfy 

Re(^2p(A)) < ... < Re(zQ 4 (A)) ^ Re(zQ 3 (A)) ^ Re(zQ 2 (A)) ^ Re(zQi(A)), (2.7) 

Re(zO,(A) - zfl j+2 (X)) > a\z\, A ^ 0, all 3 G N 2p _ 2 , (2.8) 

where 

a = 2c„_i sin — > 0, 2 = A 1/2p G 5 = (2 G C : arg 2 G — , — IV (2.9) 
4p I V zp 2pi ) 

Proof. Assume that (12.71) . (12. 8p hold for ImA ^ 0. Then identities (12. 5p give these estimates 
for Im A < 0. 

We will prove (|277]) . (jUD for ImA ^ 0, i.e. < argz < ^. Identities (Q yield 

Re(zfi p+ j(A) — 2;fi p+ j + i(A)) = Rez(a;p + j — w p+J+ i) = 2cj(—l) j Im(zrjj) (2-10) 
for all j = —p + 1, — p + 2, p — 1. Identities 

sin arg z , if 7 is even 
(-iy lm(zv 3 ) = \z\ { S ' ^0 (2.11) 

J l sm ( 2 argz) , 11 j is odd 

yield estimates (12. 7p . Furthermore, identities ( I2.10p and estimates (jl.lfip imply 

Re(zQ p+j (X) - zVt p+j+2 {X)) = Rez(u p+j - u p+j+1 ) + Rez(u p+j+ i - u p+j+2 ) 

= (-iy2 Ci Im(^) + (-l)^ +1 2c i+1 Im(^ i+1 ) ^ 2c J+1 ((-iy Im(z^) + (-iy +l lm(zrj j+1 )). 
Identity ( 12.1 ip and estimates (ll.ISp imply 

Re(zfij(A) — zQj +2 (X)) > 2c p _i max jsinargz, sin^- argzj ||z|, 

which yields (12.81) . ■ 

We define the branches rj 1 , j G N 2p , of the multiplier function r° = e tXl/2? for the unperturbed 
operator H° in the upper half-plane by the identities 

r 0(A) = e^, all AGC+. 
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Si 






1 / 


-27TJ, 
















5 3 " 




5 2 + X 




FIGURE 3. The C-plane for the cases p = 2 and p = 3, where £ = iu; 2 ? , w G Each 
sector is the image of the C± half-plane on the j-th sheet of the surface Jz?°. The values 
of the function t°(w) — at the points, connected by the lines with arrows, are equal one 
with other. 



For each j G N p the functions t\ 
2>p + j, where 



o 

p+P 'p-j+ 



T p-i+i are single-valued analytic functions in the domain 



y p+3 



C\ 



all j G N p . 



C \ M_ for even p 
C \ M+ for odd p 



(2.12) 



Lemma 2.2. i) The unperturbed multipliers satisfy the identities 



rf (A) = Tf(A), r°(A) 



r p % +1 (A) = (r^-CA))- 1 , 
Lei r° fc (A) = r° + (A) /or some A G <i 



a// (j,A) e N 2p x (C\R), 
\ R, ^ k < j ^ p. T/ien j 



(2.13) 

(2.14) 
A; + 1 and 



-MO or A = r^ n — iO for some n£N. Moreover, 



:r p+fc+ll r fc,n" 



-ZO): 



:7 p+fc+ll r 





k,n 



■iO), 



k g n;_ x 



{0,l,...,p-l}. 

(2.15) 

Proof, i) The Riemann surface Jzf of the function r° coincides with the Riemann surface 
of the function \ 1 l 2p with the unique branch point at A = 0. This surface has 2p sheets 
^jij G f%p, corresponding to the branches t®, that is Tj(X) = T°(w)\ we ^o, where A G C is 

the projection of the point w G Jz?°. For each j G N p the projection of the sheet ^ p+ j (and 
Jz^°_ J+1 ) is the domain @ p+ j and has the cut along the real axis or semi-axis, see (12.121) . The 
upper (lower) edge of each cut on the sheet Jz^° +J -, j G N p , is attached to the lower (upper) 
edge of the corresponding cut on the sheet Jz^° +J - +1 . Similarly the upper (lower) edge of each 
cut on the sheet Sf p °_j,j G N p , is attached to the lower (upper) edge of the corresponding cut 
on the sheet Jf p _j +1 . 
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The simple parametrization of the surface Jzf is given by the analytical mapping W having 
— 

the form w — >■ ( = iw 2 p , where w G ££ . We have 

W{&°) = C. 

Describe this parametrization in more details. Introduce the sectors (see Fig. [3]) 
S+ = { 2 GC:0<argz<|-}, S+ = u> j S + , Sr = {( e C : C G j G N 2p . (2.16) 

Let Jz^ 0,± , j G N 2p , be the half-plane C± on the sheet Jz^° of the surface J$f°. Then 

^(^ J ' ± ) = S±, W^ ) = all j G N 2p , where ^ = S+ U Sr. (2.17) 

Each function r°,j G N 2p , satisfies the identity t°(A) = T°(w)\ we #o, where A G C is the 

J J j 

projection of the point w G Jz?°. Let ( = W(w) G C,w G Then r°(w) = e*>. Using 

identities (I2.17P we obtain r?(A) = e^, where is given by the conditions Q G £j,A = 

(-l)% 2p . Then A = (-1) P C/ P (see (12~T6|) and Fig. E} and the identities r°(A) = = = 
r°(A) for all (j, A) G N 2p x (C \ R) give the first identity in f[2TT3]) . Identities (J23H imply the 
second identity in (gUS]) . Identities ([USD Y ield (EUD- 
ii) Let r p ° +fc (A) = r p ° +i (A) for some AeC\R,0<A:<j<p. Then 

e C P+fe=e C P+J; C 2^ = ^ = ( _ 1)PA> Cp+fcG5p+fe; (p+j . eSp+ , (2.18) 

The first identities in (I2.18P yield ( p +k — ( p +j = i^nri for some n = ±1, ±2, ... Then we have 
Re( p+ k = Re(p+j- The second identities in (12.181) imply \( p +k\ = \Cp+j\ = I z l> which yields 
Cp+fc = Cp+j- Moreover, lm( p+ k — lm( p+ j = Inn. Identities (12.161) and the condition k < j 
(see also Fig. 0J) give: a) j = fc + 1; b) ( p+k = f e^, i.e. A = = (-l)*(^) 2p = r° n ; 
and c) identities (I2.15p . ■ 
Consider the operator = (— l) p Jp^ + H ^pX ,^ G R- The equation 

(-l)V 2p ) + /V 2p_2) = Ay, (t, A) G R x C, (2.19) 
has the solutions e u i^ zt ,j G N 2p , where Uj(X),j G N 2p , are the solutions of the equation 

( W M)2 P + (_l)P( w M)2 P -2 /i ^2 _ = Q (2>20) 

The functions u)j,j G N 2p , constitute branches of the analytic function uj^ 1 having only alge- 
braic singularities. For A G Ar for some R > large enough, we define these branches by the 
asymptotics cuf (A) = + o(l) as |A| — > oo, here and below 

A R . = {A G C : |A| > i? 2p }, i? > 0, A± = A fi n C±. 

The multipliers have the form t^(A) = e 2 ^ ^ for all (j, A) G N 2p x for some R > large 
enough. Then 

7f(A) = e^W, all (j,A)GN 2p xA fl , where n?(A) = H^' ° ,(2.21) 

I cj^ A J , 1m A < U 

where £ is given by identities (12. 6p . Each function (A), j G N 2p , /i G R, is analytic in A G A fl 
and piecewise-continuous in A G C, but the set {f^(A), j G N 2p },/i G R, is continuous in A G C. 
The branches of the Lyapunov function are given by A^(A) = coshzf^A), (j, A) G N p x A R . 
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Lemma 2.3. i) Each function u^,j G N 2p; satisfies the asymptotics 



u)f(\) = Uj - )- — + 0(|2:|~ 4 ) as |A| -> oo, (2.22) 

J 2pL0jZ z 



^(A + e) = ^(A) + 0(|^|- 4 ) as |A| ->• oo, e = 0(|^| 2 ^ 2 ). (2.23) 
TTie periodic and antiperiodic eigenvalues for equation (— l) p ?/ 2p ) + /i?/ 2p ~ 2 ) = Ay satisfy: 
A£' + = 0, A£~ = A£'+ = (vm) 2p - (-l) p ^(vrn) 2p - 2 , all n > 1. (2.24) 
raj £^ac/i function G N2 P , satisfies the asymptotics 

Q»(\) = Q J (\) + 0{\z\- 2 ) as |A|^oo, (2.25) 

Proof, i) Substituting ojf = uj + 5, 5 = o(l), into the identity (t^) 2p + (-l) p (u^) 2p ~ 2 fiz' 2 - 
(— l) p = we obtain 

{Uj + 5) 2p + (-l) p ^j + 5) 2p ~ 2 fiz- 2 = (-l) p . (2.26) 

This identity gives 2puf~ 1 5 + 0(5 2 ) + 0(\z\- 2 ) = 0, which yields 5 = 0(\z\- 2 ). Using 
(I2.26P again we obtain 2pu 2p ~ 1 5 + (— 1) p uj 2p ~ 2 fiz~ 2 = 0(|z|~ 4 ). This asymptotics gives 5 = 
(-l)p+ l (2pu j )- 1 fiz- 2 + 0{\z\- A ), which yields (I2T22]) . Asymptotics fl2T22|) yield 



^(A + e)-<(A) = ^-^(---) +0(|^|- 4 ) as |A| ^ oo, (2.27) 



where C = (A + e)2F. Using the asymptotics £ = 2 + OQzl -1 ), we obtain (I2.23p . 

ii) We will prove ( I2.24p for the periodic eigenvalues. The proof for the antiperiodic eigen- 
values is similar. The periodic eigenvalues Aq' + , A^, n ^ 1, are zeros of the entire function 
Dt = n?=i(Ai - !)■ We have 

P P M \ P 00 / I An2 2\ 2 

o?=n( c ^-i)=2'n rf ^=^riK) 2 n( i +7&) ■ ^ 

j=l j=l j=l n=l V V I ' 

Using the simple identity 

p 

u 2p + {-\) p uj 2p - 2 iiz- 2 - (-l) p = Y[(u 2 - (c^) 2 ), all wGC, 
we obtain n^iO^) 2 = ~~ 1 (P u ^ w = 0) and 

n(^)-^n((^)'-^)- ( - * 

Substituting these identities into (12.281) we obtain 

A frA (-1)> A 



(2mr) 2 (2n7r) 



2p" 



^"LiV (2mr) 2 (2mr) 2p ^ 

which yields (I2.24p for the periodic eigenvalues. 

iii) Asymptotics ( jX22D and definition (?2~2~T\> of flj yield flM 
Remark. The periodic eigenvalue A = for equation (I2.19P is simple and other periodic and 
antiperiodic eigenvalues have multiplicities 2. 
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3. Fundamental solutions 



In this section we consider the operator H = H + q, H = (— l) p Jp^ and recall that 
? = qj e Ll eal (T), jeN p = {l,2,...,p}, p^2, 

3=0 

The form domain of the self-adjoint operator Hq is the set Domj-^o) = W 2 (WL). The 

quadratic form (qy,y) is defined by (qy,y) = Y^=o(~ ^YiQj+W^i V^)i V e Wp(R). Let 

1j = Qj,o + a p where g^o = qj{t)dt and aj G Wi(T), J a,j(t)dt = 0, j G N p . The integration 
by parts gives the form 

p-i 



(qy,y) = U^o\\y U) \\ 2 - 2Re(a,W i+1) ,2/ (i) )) 

3=0 ^ ' 



(3.1) 



correctly defined on the form domain y G W 2 (R). 

For each , j G N p , we introduce the sequence (qj, n )%Li of the smooth functions q^ n G 
VI^L^T), such that 

/ {o.j ~~ Qj,n)dt = all n G N, /3 n = sup / |gj — qj, n \dt — > as n — >■ oo. (3.2) 
Jo ieNp Jo 



Let 



^ d j d j 

H n = H + q {n) , neN, where g (n) = ^ l+i q ^ n l+r ( 3 ' 3 ) 

j=o 

Proposition 3.1. The quadratic form (qy,y) satisfies 

l(^,y)l^^l|y (p) H 2 + C||y|| 2 , all yeW 2 (R) (3.4) 

/or some constant C > 0, where \\y\\ 2 = (y,y) is the scalar product in L 2 {R). 

ii) There exists a unique self-adjoint operator H = H +q with the form domain Dom^(iJ) = 
W 2 (R) and 

(Hy,y 1 ) = (H y,y 1 ) + (qy,y 1 ), all y, y x G W 2 (R). (3.5) 

Hi) Let y G W 2 (M.). Then there exists the sequence (e n )f > such that e n > for all neN, 
£ n — > as n — > oo, and 

\((H-H n )y,y)\^e n (\(Hy,y)\+\\y\\ 2 ), all neN. (3.6) 

Proof, i) We get 

2\{a j y^\y^)\^2b\\y^\\\\y^\\ < e\\y^f + -\\y^\\ 2 , b= sup \ aj {t)\ (3.7) 

£ (tj)6TxN p 



for any e > 0, where (y, y) = \\y\\ 2 = L \y\ 2 dt. Thus (13. ip yields 

\(qy,y)\ < £ ||^|| 2 + 5^ ||yW|| 2 , g = q + e + % = ^ |g i)0 |- (3.8) 
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Using the simple estimate Ylo 1 lls/^ll 2 ^ £2 ||z/^l| 2 + Cer 1 1 2/ 1 1 2 f° r some C £ > 0, we obtain 

\(qy,y)\^e(l + eC)\\y^\\ 2 + CC E \\y\\ 2 , (3-9) 

which yields (j3.4p with | = e(l + sC) and C = CC £ for e small enough. 

ii) The operator q on the domain Dom.f d (H ) is given by (13.11) and satisfies the estimate 
( 13. 4p . Using the KLMN theorem (see |RSlj ) we obtain that there exists a unique self-adjoint 
operator H = H + q with the form domain W 2 (M.) and identity (13. 5p holds true. 

iii) Repeating the previous arguments and using ||?/ jf ' ) || 2 ^ ||y || 2 + \\y\\ 2 we obtain 

^ 2/3 n (||^)|| 2 + \\y\\ 2 ) = 2f3 n ((H y,y) + |M| 2 ), 
where (3 n is given by ( 13. 2p . This yields 

p-i 

\({H-H n )y,y)\ <: ^ |((g i+1 - q j+1>n )y^ ,y^)\ <: 2p(3 n ((H y,y) + \\y\\ 2 ). (3.10) 
o 

Using estimates (13.41) we obtain 

(H y,y) ^ \{Hy,y)\ + \(qy,y)\ ^ \(Hy,y)\ + ^(H y,y) + C\\y\\ 2 , 

which yields \(H y,y)\ ^ 2\(Hy, y)\ + 2C\\y\\ 2 . Substituting this estimate into ( I3.10p we obtain 
( BSD , where £ n = 4pf3 n {l + C). m 

Below a vector h = (h n )f G C N has the norm \h\ = \h n \, while an iV x Af matrix 
A = (Aij)f^j =1 has the operator norm given by |^4| = supiM =1 \Ah\ = maxi^j^ YuiLi l^jl- 
Always below we denote n x n diagonal matrices by diag(aj)™ =1 = (ajO~jk)^k=i- The following 
Lemma (proof in Appendix) describes the basic properties of the monodromy matrix A). 

Lemma 3.2. The matrix-valued function Ai(l,X) is entire and satisfies: 

|A4(1,A)| < 2pe zo+ *, all XeC, \Z-\\)M{1, X)Z{X)\ < 2pe Z0+K , a/Z |A| > 1, (3.11) 

where 

Z = diag(z J_1 ) 2p , x = max / \q~(t)\dt, zn = max I Re(zu) 

jen p J Q j<=N P 

Moreover, Ai(l,X) is a continuous function of (qj){ G L 1 (T) P . 

Lemma 3.3. Let each qj G Wj_j_(T), j G N p and let A G C. Then the spectrum of the 

2p x 2p matrix (ip^ (1, A))^ P J=1 coincides with the spectrum of the matrix X), counted 

with multiplicity. Moreover, in this case 

a(H) = {X G R : |r,(A)| = 1 for some j G N p }. (3.12) 
Proof. Introduce the vector- valued function Y = (y^~ 1 )) J w L . Identity ( II ,5p shows that 

/0 ... q p 0\ 

(3.13) 



Y = SY, where S = [ V I , S 



S 1 



p-i 



... 

\0 q 2 ... 0/ 
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Let M(t, A) = (ipf 1} (t, X))l p j=v (t, A) G IxC. Identity f lBTTgj) shows that each matrix-valued 
function SM(-, A), A G C, satisfies equation 0Eg, which yields M(t, A) = S(t)M(t, A)S -1 (0). 
Using S{1) = S{0) we obtain M(l, A) = S(0)M(1, A)«S _1 (0)- Thus the matrices M(l, A) and 
Ai(l,X) are similar and the spectra of these matrices coincide one with other, counted with 
multiplicity. Identity fl3~12]) follows (see, e.g., [DS], Th. XIII.7.64). ■ 

We will introduce the Jost type fundamental matrix solution T(t, A) of equation (13.171) . This 
solution will be described below. Recall that the monodromy matrix has the form -M(l, A), 
where the matrix- valued function A4(t, A) satisfies the matrix equation (II. 4p . Rewrite equation 
(ll.4p in the form 

M' - V fM {X)M = Q**(t)M, (t,A)6lxC, (3.14) 
where the 2p x 2p matrices V 1 and Q M are given by 



V» = V- (-!)>£, Q 1 = Q + (-!)*>£, n 



q p (t)dt G 



the matrices V, Q are defined by (11.61) . and the matrix £ is given by 



£ 



1, £ jk = 0, all (J, k) ^ (p + l,p). 



(3.15) 



(3.16) 



Each matrix P M (A), A G A/j = {A G C : |A| > R 2p } for some R > large enough, has eigenval- 
ues zQj(X),j G N2 P , all these eigenvalues are simple and the corresponding eigenvectors are 
given by 

/ 1 \ 



U: 



W) 2 



{zwy +l + i-iy^zny- 1 



3 en 



2p- 



\(^) 2p - 1 + (-i)>(^) 2p - 3 / 

Then the matrix V 1 is similar to the diagonal matrix 

zB" = U- l Vm, where = diag(^) 2p , 

and the 2p x 2p matrix U has the form 

w = {u x u 2 ... u 2p ) . 

We rewrite equation (13.141) in the form 

M' — zB^(X)M = Q(t, X)M, (t, A) eRx A r , 
where 

M = U- l MU, Q = {Q tj )Z=i = U- X Q?U. 

The following Lemma, proved in Appendix, shows that the matrix Q is decreasing at large 
IAI. 



(3.17) 
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Lemma 3.4. The matrix-valued function Q satisfies the following asymptotics: 

Q(t,X)= { ^^((q P (t)-^C(X)+b(t)0(\zn), C = {C jk )% =l , C jk = ^-, (3.18) 

uniformly on t G [0, 1] as |A| — > oo, where b(t) = max^p l^jCOl- 

Consider equation ( I3.17p . Assume that this equation has the 2p x 2p matrix- valued solution 
T(t, A) for all t G [0, 1] and some A G A R . Then M(t, A) = T(t, A)T(0, A)" 1 and 

M(t, A) = U(X)M(t, \)U-\\) = U(\)T(t, A)T(0, A)- 1 W~ 1 (A), (t, A) G [0, 1] x A R . (3.19) 

In order to analyze equation ( 13. 17[) by the Birkhoff method (see [Na] ) . we write the solution 
T in the form 

T = Ge zBH , (3.20) 

where Q(t, X),(t, A) G [0,1] x A R , is some matrix-valued function. Substituting (I3.20p into 
(I3.17P we obtain 

G' + z{GB» - B^G) = QG. (3.21) 
This equation is equivalent to the integral equation 

G = l 2p + KQG, (3.22) 
where K is an integral operator given by 

f Joe z{t - s){n ^ x] -^ W) A ij (s)ds, i£i>j f 1 
(KA)ij(t, A) = < _ °i e ^(t- s )(nr(A)-n^(A))^^ s ^^ if i < j = i ^ ~ S ' A )^( s ) rfs ( 3 ' 23 ) 

for the matrix- valued function ^4, where (i,j,t, A) G N| p x [0, 1] x A R and 

A ) = I^w^k if i < ' x(t) = (i, if t > o • (3 - 24) 

In fact, differentiating fl 3 . 2 2 j) we obtain 

= (QG) ij (t) + z(^(X)-^(X))(KQG) ij (t 1 X), all (i,j,t, A) G N|p x [0, 1] x A R . (3.25) 
Using the identities 

(Of - = (B^KQQ) - (KQG)B% = {B»G - GB% d 

we obtain G' = QG + z(B^G — GB^). Thus G satisfies (13.211) . which yields the equivalence of 
equations fl3~2TD and (137221) . 

In Lemma \'3. 5 1 proved in Appendix, we describe the Jost type fundamental matrix solution 
T of equation (13.171) . In fact we will prove that equation (13 .22 j) has a unique solution for |A| 
large enough. For the 2p x 2p matrix- valued function A G L°°(0, 1) we introduce the norm 

, = sup \A(t)\. 

te[o,i] 
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Lemma 3.5. i) Let A G L°°(0, 1) be a 2p x 2p matrix-valued function. Then for each A G Ar 
for some R > large enough the operator K satisfies the estimate 

WiKQA^-.X)^^ ||| where £ = ™ax{4^ - fi\dt; l}. (3.26) 

For eac/i A G £/ie integral equation §3.22$ has the unique solution Q(t,X). Each 
matrix-valued function Q(t, -),t G [0, 1], is analytic in A^ and satisfies the estimates 

\\g(-,X)\\ 00 ^2, \\g(;X)-t 2p \\oo< %, ||S(-,A)-l2p-0i(t,A)||oo<S, (3.27) 

/or all X G A 2 #, where Q\ = KQ, 

(Gi(t,X)) ij = ' 3 \ e i3 (t-s,X)(q p (s)-fi)ds + 0(\z\- 2 ), i,j G N 2p , (3.28) 

Jo 

as \X\ —> oo, uniformly on t G [0, 1]. 
Hi) For each t G [0, 1] t/ie function T(t, •) zs analytic in A^ and satisfies 

T{t,X) = e zBfl W t {l 2p + 0{\z\- 1 )) as |A| -> oo, (3.29) 

uniformly on t G [0, 1]. 

Now we will prove the main result of this Section. In this Lemma 13.61 we obtain the repre- 
sentation and asymptotics of the monodromy matrix A4, see fl3.30l) - (l3.34l) . 

Lemma 3.6. i) The monodromy matrix •) satisfies the identity 

M{l,-)=Ug{0,-)Fe" B,t (Ug(0,-))~ 1 , where 7 = Q-\0, -)£(1, •)• (3.30) 

ii) The matrix-valued function J 7 = (J-ij)^ =1 is analytic in A^ and satisfies the asymptotics 

^i(A) = % + 1 - ; Hfof 1 / A)( 9p (t) - + 0(|^|- 2 ), (3.31) 

.Fy(A) = Od^r 1 ) /or i ? j, JT. .(A) = 1 + 0(|^|- 2 ) (3.32) 
as |A| — > oo, where 

r e 2(i-t)(n i (A)-^-(A)) ) ^ > j 

Hi) The functions F p+ k, P +k+i and Fp+k+i^+k for all k G N?_j = {0, 1, ...,p — 1} satisfy 

J r p+feiP+fc+ i(A) = / fc)n + 0(n~ 2 ), J r p+fe+ i >p+fc (A) = 7 fci „ + 0(n~ 2 ) (3.34) 
as n -> oo, A = (-lW™)2 P + 0(n 2p " 2 ), Im A ^ 0, w/iere 

2(-l) fc+1 c fc J e^g p>n , fc odd 
2p7rn (^-e 2 fg p n , fc even 

Proof, i) Identities flQ5|) . flg^Oj) yield M(*,A) = U(X)G(t, X)e zBn g- 1 {0, X)U-\X), which 
implies (I3.30p . 

ii) Estimates (I3.27P yield {?(t, A) = 1 2?3 + A) + 0(|2:|~ 2 ) as |A| — > oo, uniformly on 
t G [0, 1], and 

F(\)=g- 1 (0,\)g(l,\) = l 2p + g 1 (l,\)-g 1 (0,\) + O{\z\- 2 ) as |A|->oo. 
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Substituting ( I3.28P into this asymptotics we obtain 

Jij(A) = Sij + / ( eii (l - 1, A) - eo-(-t, A))(g p (t) - u)dt + 0{\z\- 2 ) 

^P z Jo 

as |A| — > oo. Substituting ( I3.24p into the last asymptotics and using (I2.25P we obtain f )3.3ip . 
which yields the first asymptotics in (I3.32p . The second asymptotics in (I3.32p follows from 
the identities ^ = 1, see (13. 33 p . and the identity q p (t)dt = n, see (I3.15p . 

iii) Identity gives = w i for all j G N 2p . Identities Q, £fc 1 = £fc and ef = (-l)P+ fc 
(see (T23D) yield for jfe G Nj-i= 

ulul-l = = (-l) s e kVk , uf s+1 ur l = e 2 r l Vk = (-l) s e kVk , s = p + k. (3.36) 

Substituting identities (I3.36P into asymptotics (13.31 j) and using (I3.33|) we obtain 

^ >s+ i(A) = { - 1)k+1 " kT]k f e ^-^\q p (t) - rfdt + 0(\z\-% 
A P Z Jo 

•^+1 >S (A) = { ~ l)k+l£kllk f e^-^-«°\q p (t) - n)dt + 0(\z\- 2 ) (3.37) 
^P z Jo 

as |A| -> oo, G Nj-i- Let A = (-l) fc (^) 2 ? + 0(n 2p ~ 2 ) as n ^ oo. Then z = A^ = 
2^ + 0{n- 1 ) and using (Q we have z(u s - w s+ i) = i(-l) k+1 2nn + 0(n -1 ). Substituting 
this asymptotics into ( I3.37P we obtain 

^, s+ l(A) = (~ 1 ) fc+lgfcCfc ^ e <(-l>**n* ((fc(f) _ ^ + 0( „- 2)j 

^ + i, s (A) = Ljj!!^ ^-D^ 1 *™*^) - ^ + 0(n- 2 ) (3.38) 
2p?m y 

as n — > oo. Substituting e k from (I2.4[) into (I3.38|) . we get (13.34j> . ■ 

4. Properties of the multipliers 

Define the single-valued branches of the multiplier r at high energy. Here we use the results 
of Lemma I4T31 which will be proved later. The zeros of the function p for high energy are close 
to the real axis. Then the functions Tj,j G N2 P , are analytic in the domain Ar R {A G C : 5 < 
arg A < 7r — 5} for some R > large enough and for any 5 > small enough. Asymptotics 
( II .8p define the branches Tj,j G N2 P , of the function r in this domain. 

Moreover, the points r^ n G A R , k G {j — l,j},j G N p _i,n ^ n for some (large) n ^ 1, 
are ramification points of the functions t p+ j and T P -j+i and these functions have no any other 
singularities in Ar (see discussion after the proof of Lemma [4.41) . Here r^ n satisfy: 

!) \ r tn- r ln\ <1 f o r all k G N°_ x = {0, ...,p — l},n ^ n , 

2) all TQ n ,n ^ no, are positive numbers (anti-periodic and periodic eigenvalues), 

3) allr± n ,fc G N p ,n > no, are real or non-real numbers (ramification points of the Lyapunov 
function), ± Imr^ n ^ and if Imr^ n > for some k G N p , n ^ n , then r^ n = f^ n . 

Then each function t p+j -, r p _ :)+1 , j G N p , is a single- valued analytic function in the domain 
@ p +j, where 

% +j = a r \ u n>no (Tj. ltn u r jt n) ail j g Np_i, % p = a r \ u n>no r p „ hn , 
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each Tj )Tl jj G Np_ 1; n ^ n , is the segment \r^ n _ x , rj >n ] of the line. 

Note that for each j G N p the cuts r° n = [ry n _ x , ], n ^ 1, on the real axis for the 
multiplier function r° of the operator H° merge and constitute the cuts along the negative 
(for odd j) or positive (for even j) semi-axis, see (I2.12p . 

The following result is a simple consequence of the standard perturbation theory. 

Lemma 4.1. i) Let X G for some R > large enough. Then for each j G N 2p the 
monodromy matrix A) has the eigenvalue Tj satisfying the asymptotics 

Tj (\) = e zn >(l + 0(\z\- 1 )) as |A|-»oo, (4.1) 

which yields M.8\) . Moreover, if Tj(X) = r k (X) for some 1 ^ j < k ^ 2p and some X G A R , 
then k = j + 1. No more than two multipliers Tj,j G N 2p; can coincide one with other at the 
point X. 

ii) The multipliers satisfy the identities 

r P-i+i(A) = Tp+jW, all (j,X)eR p x% +r (4.2) 



Proof, i) Recall the following Gershgorin's result from the matrix theory, sec [HJJ: 

Let A = (v4jj)™ =1 be a complex m x m matrix and let Ri = Y^JLij^i f or * £ N m . Then 
every eigenvalue of A lies within at least one of the discs {r G C : |r — Au\ < Ri},i G N m . // 
the union of k discs is disjoint from the union of the other n — k discs then the former union 
contains exactly k and the latter n — k eigenvalues of A. 

Identity (13.301) implies that the matrices -M(l, ■) and X = e z T have the same eigenvalues 
(the multipliers). Asymptotics (I2.25p . (13.321) give 

X = e zB ° ^l 2p + ' where B ° = di ag(%)f , 

and the matrix-valued function W(z) = {wij(z))fj =1 is uniformly bounded on \z\ > R for 
some R > 0. By Gershgorin's theorem every multiplier r lies within at least one of the discs 

Xj = |r G C : |r - e zQj \ < \e znj \-^ j, j G N 2p , where wj = max^ \w jk (z)\. 



Firstly, let the disc ^ for some j G N 2p be disjoint from the other discs Jlf k) k ^ j. Then, by 
Gershgorin's theorem, the disc Jtfj contains exactly one multiplier tj(z), which satisfies the 
estimate \tj(z) — e zQ: >\ < This estimate gives (14.11) for this case. 

Secondly, consider all k, j, 1 ^ j < k ^ 2p, such that J^j fl J^ k ^ 0. Then the distance 
between the centers of these discs is less than the sum of their radii: |e 2flj — e zUk \ < (|e 2 ^ \wj+ 

\e zUk \w k ). Then 

le z(n k -n j} _ X | < 1/ + (4 _ 3) 



where we used the estimate | e z ( nfe -^)| = e Re«(n fc -^) ^ ( gee (j^). U k ^ j + 2, then 
estimates (|2T8|) together with (LTTj) yields Re z(Q k — Qj) < -a\z\,a > 0, and then e z ( n k-n 3 ) q 
as \z\ — )■ oo. For |z| large enough we have a contradiction with (14.31) . Hence k = j + 1. 
Moreover, the similar arguments show that only two domains J^, J?^+i can intersect each 
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other and they are disjoint from other domains J^, m ^ j,j + 1. In fact, let Jfcj PI Jfrj+i ^ 
and Jtj + \ fl J£j+2 7^ 0- Then we have two estimates 

|e*nj_ e *ni+i| < n e ^\ Wj+ \ e ^+i\ Wj+1 ^ | e ^+i_ e ^i+2| < ^(| e ^+i| u , j+1 + | e ^+2| lt ; i+2 ) ) 
which yield 

\e zn i -e zn i+ 2 \ < | e ^-e z ^ +1 | + |e^ +1 -e 2 ^ +2 | < ^{\e z ^\w j + 2\e z ^ +1 \w j+1 + \e zQ]+2 \w j+2 ) ■ 
Then 

|1 _ e *(n,-+a-«i)| < J_( Wj + 2Wi+1 + Wj+2 ), 

which is in contradiction with the estimate Re z(Qj +2 — < — a|z|,a > 0. Thus only 
two domains Jfrj, J^+i can intersect each other and they are disjoint from other domains 
Xn,m=£ j,j + 1. 

Let J£j PI 7^ for some j G N 2p _i. By Gershgorin's theorem, the domain U J^ + i 

contains exactly two multipliers Tj(z), Tj + x(z) and |t,- +1 — e 2 ^ +1 | < -^(2\e znj \wj + \e znj+1 \wj +1 ), 
which implies 

\r j+1 e- z ^ +1 - 1| < |3(2tWj + |e z(n '~^ +l) |wj+i)- 

Estimate (14.31) yields that e ^(^-^j+i) = l + Od^l -1 ) and then r j+l satisfies asymptotics (14. ip . 
The similar arguments show that r,- also satisfies (14. ip . In fact, jr^ — e 2 ^| < jj|(|e zf ^|u>j + 
2|e which implies 

| Tj e"^ - 1| < ^(2wj+i + le^'+^^Vi+i)- 

Since e Rez ( n j+i-^j) ^ 1 ; we obtain asymptotics (14.11) for Tj. 

ii) Asymptotics (14.11) and identities ( I2.14p provide ()4.2p . ■ 
Proof of Theorem 11.11 The proof of identity (ll.lOp is standard. We rewrite equation (jl.4p 
in the form J M! = 'KM., where T-L = J(V + Q), JQ is a diagonal 2p x 2p matrix and 

Then the matrix "H is symmetric: "H T = H. Using the identity J~ T = —J we obtain 
— {M^yj = M. J, H, which yields 

(M T JM)' = (M T )'JM + M T JM' = -M T UM + M T UM = 0. 

Then M T JM = const and using M(0) = l 2p we obtain (ll.lOp . 
i) We will use the arguments from |CK] . Identity ( ll.lOp yields 

D(r,-) = r 2 ^(r- 1 ,-), r^O, (4.4) 

and -D(r, A) = Ylo P > c k{^) T ' 2p ~ k i where the functions Xk are given by 

1 2 *Vt_lt ^ Tr.M fc (l,-) 
x = 1, xi = -2pT u x 2 = — -(T 2 +TiXi), x fc = — — y Tk-jXj, T fc = — , 




o y 



EVEN ORDER PERIODIC OPERATORS ON THE REAL LINE 



21 



(sec [RS2J, p. 331-333). By Lemma I3T2") the coefficients x fc (A) are entire in A G C. Using the 
identity (@~1D we obtain x 2p -j = G N°, which yields D(r, •) = (r 2p + 1) + xi(r 2p_1 + r) + 
... + x p _i(r p+1 + r^ 1 ) + x p r p . Then 

= ^ + /i(A)f p_1 + ... + / P (A), ,/ = I±^, (4.5) 

where /i, .., f p are some linear combinations of x , .., x p . In particular, all coefficients /i(A), 
/ P (A) are entire functions. The function $(z/, A) = is a polynomial of z/ of degree p. Each 

zero Aj, j G N p , of this function satisfies Aj = \{Tj+rJ l ), j G N p , where Tj, tJ are multipliers, 
and identity (11. lip holds. Asymptotics ( 11.8)) yields ( I1.12p . Recall that Tj are branches of the 
function r analytic on the 2p sheeted Riemann surface (see Sect.l). Then Aj(\),j G N p 
constitute p branches of one analytic function A(A) on the connected p-sheeted Riemann 
surface ffl. 

ii) Proof repeats the standard arguments (see [CK] ). 

iii) Let Ai n ,n ^ 1, be the monodromy matrix for the operator H n given by (13.31) . Let 
Tj >n , j G N p , be the multipliers of H n . Identity (I3.12p gives 

a{H n ) = {A G R : |r i)n (A)| = 1 for some j G N p }. (4.6) 

Lemma I3T21 provides M. n — y M.(l, ■) asn — y oo uniformly on any compact in C. Then Tj >n — > Tj 
as n — > oo uniformly on any compact in C for all j G N p . Identity (14.61) implies 

a(H n ) -)> {A G R : |r,(A)| = 1 for some j G N p } as n -)> oo. (4.7) 

Assume that 

er(# n ) -> er(if) as n -> oo. (4.8) 
Then using (14 .7p we obtain cr(-ff) = {A G R : |Tj(A)| = 1 for some j G N p }, which yields (11. 13j) . 
Now we will prove (14. 8p . We need the following result (see [Kaj . Th. VI. 5. 13, Cor. V.4.2): 

Let be an operator in a Hilbert space Jtf? and let Q n be a symmetric operator in ffi 1 
with the form domain Dom.fdQ n C Dom^ A. Assume that 

\(Q n u,u)\ ^ e n (\\u\\ 2 + (Au,u)), ueDom fd Q n , where e n > 0, e n -> 0. 

Then the the Friedrichs extension A n of the operator A + Q n is self adjoint for sufficiently large 
n and A n — > A in the uniform resolvent sense. If a (A) has a gap at a, then cr(A n ) has a gap 
at a for sufficiently large n. 

Estimate (13.61) shows that H n — y H in the uniform resolvent sense (and then in the strong 
resolvent sense) as n — > oo, and lim^oo cr(H n ) C cr(H). Then relation (14. 8 p is obtained from 
the following result (see [Ka], Th. VILLI. 14): 

Let H,H n ,n ^ \, be self adjoint operators in a Hilbert space ffi and let H n —y H in the 
strong resolvent sense. Then every open set containing a point of cr(H) contains at least a 
point of a(H n ) for sufficiently large n. ■ 

Consider two simple examples. 
Example 1. Consider the operator H = (— 1) P J^ + X]j=o<?j+ijpJ with the constant coeffi- 
cients <2j. Equation ( I1.3P has the solutions e ±? ^^*, where £j(^) = \f w j{^)i3 e N p , and Wj 
are values of the algebraic function w(X) which is a solution of the equations P(w) — A = 0, 
P(w) = w p + Y^=o(~ Qj+i w '' ■ Then the multipliers have the form and the Lyapunov 

function is given by Aj(A) = cos£j(A). If the polynomial P(w) — A for some A G R has n 
(1 ^ n ^ p) positive simple zeros, then the spectrum a(H) in some interval (A— s, X+e),e > 0, 
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has multiplicity 2n. Let P(w) = T p (w - 1), where T p (w) = E? /21 {L)( w2 ~ l) n w p ~ 2n is the 
Chebyshev polynomial. The properties of these polynomials (see |ASJ) provide that the spec- 
trum is given by cr(H) = [— 1, +00) and has multiplicity 2p (maximal multiplicity) on the 
interval (—1,1) and the multiplicity 2 (minimal multiplicity) on (1, +00). 
Example 2. Consider the operator H = H p , where H = —4$ + q is the Hill operator 
with the 1-periodic function q G L 2 (T). The branches of the Lyapunov function are given 
by Aj(A) = A(—z 2 u 2 ),j G N p , where A is the (entire) Lyapunov function of H. Recall 
that the spectrum cr(H) is semi-bounded below and consists of bands separated by gaps. Let 
<jj = {A G E : A(-z 2 uf) G [-1, 1]}, j 6 N p ,z = A* G S. The spectrum of H satisfies the 
identity cr(H) = U^o"j = a± U a p . We have 

{{A G [0, 00) : —z 2 G cT-i, p even r , r ^ . 9 . ~ r ^ 

f ,\ 1 ' ' \ , <7 p = {Ae 0,oo :* 2 e<r(# n 0,oo }, 

{A G (— 00, 0J : — e p z G cr_}, p odd 

where <r_ = cr(if) fl (— 00, 0]. If p is odd, then the spectrum <r(H) has multiplicity 2. If p is 
even, then the spectrum has multiplicity 4 in the set <j\ fl cr p and multiplicity 2 in the other 
intervals. 

Recall that all functions p,D± = 2~ P D(±1, ■) = ]T^ =1 (Aj =F 1) are entire. The zeros of 
p are ramifications of the Lyapunov function, the zeros of D± are periodic and antiperiodic 
eigenvalues. We introduce the contours C n {r) = {A : \z — irn\ = rrr}, r > 0, n ^ 0. 

Lemma 4.2. Let N G N be large enough. Then the function D + (and D_) has exactly 
2N + 1 (and2N) zeros in the domain {\z\ < 27r(A^+|)} (and in {\z\ < 2nN}), counted with 
multiplicity. Moreover, for each n > N the function D + (and D_) has exactly two zeros in 
the disk {\z — 2irn\ < f } (and in {\z — 11 (2n + 1)| < \}), counted with multiplicity. There are 
no other zeros. 

Proof. We consider the function D + . The proof for D_ is similar. The function D + for the 
operator H° is given by 

i*=n(co*«*-i)=-£n(i-^)* 

j=l n=l v \ I / 

Assume that for each j G N p and for some R > large enough 



coshzco,- — 1| 



sinh 



2 >l e l Re ^l, all AGA i? \U neN {|^-27rn| (4.9) 



2 

Then asymptotics (I1.12p and estimates ( 14. 9 p yield 

£±W _ TT A 3 -(A)-1 _ fr cosh^-l + Od^l-V^^') _ , 

m\) "J- J- cosh "l-l- coshs^-l + J l^.lUj 

+v ; i=1 3 j=1 

as A G C \ U„ eN {|2: - 27rn| ^ f }, |A| -)■ 00. Let AT ^ 1 be large enough and let N' > N 
be another integer. Let A belong to the contours Co(2A" + ~),Cq(2N' + |),C2 n (|), \n\ > N. 
Asymptotics (I4.10p yields 

D + (X) 



D + (X) - Dl(X)\ = \Dl(X)\ 



D° + (X) 



- 1 



\Dl(X)\0(\zn < \D° + (X)\ 
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FIGURE 4. The domains %k,n m the z-plane for p = 6. 



on all contours. Hence, by Rouche's theorem, D + has as many zeros, as in each of the 
bounded domains and the remaining unbounded domain. Since has exactly one simple 
zero at A = and exactly one zero of multiplicity two at (2irn) 2p ,n ^ 1, and since N' > N 
can be chosen arbitrarily large, the statement for D + follows. 

We will prove ( 14. 9p . Using the simple estimate e' Im ^ < 4| sinz| as \z — Trn\ > f for all n G Z 
(see |PT] . Lemma 2.1) we obtain that estimates (14. 9 p hold in the domain C \ U ne %{\zcUj + 
i2im\ ^ §}• For each j ^ p the estimates \zuj + i2-nn\ > | hold for all n G Z and |A| large 
enough. Moreover, the estimates \zu p + i2irn\ = \z — 2im\ > - hold for all n ^ and |A| large 
enough. Thus, estimates (14. 9p hold in A R \ U ne ^{\z — 2nn\ ^ |}. ■ 

Now we will describe the zeros of the function p. Identifying the sides of the sector S = 
{z G C : axgz G (— ^,^]} (i-e. we identify each point xe l ^p,x G R+, on the ^-plane with 

the point xe~ l ^p) we obtain the cone S con . For each (k,n) G N p _i x N°,N° = N U {0}, we 
introduce the domain n given by 



k.n 



A 



t 2p ,zeS c , 



irn- 



) 



< p\, f3 > is small enough. (4.11) 



Each domain is a neighborhood of the zero r° n = (irn^) 2p of the function p°, see (I1.16p . 
where rj^ is given by (I2.4p . 

We have ty lfl = ... = %-i, and <& 1)0 n = for all (k,n) G x N. If p = 2 or 3, 
then the domains ^k,n,^k',n', (k,n), {k',n') G x N, (k',n') ^ (k,n), are separated, that 
is ^k,n H ^fc',n' = 0. The situation is more complicated for p ^ 4 (see Fig. HJ). In this case 
we have ^2fc,n H %2k'-i,n' = for all 2k, 2k' — 1 G N p _i,n, n' G N. However, the domains 

N p -i,n,n' G N such that 



%2k,ni ^2k',n' can have non-empty intersection for some 2k, 2k! G 
(k',n') ^ (k,n). The similar statement for the domains ^2fc-i,n> ^W-i.n' holds. 

We introduce the cluster decomposition of the set of indices (&,n) G N p _i x N°, having the 
form U|t_ 00 ^ = N p _i x N°, where 

TTie indices (k,n),(£,m) belong to the cluster ffj iff the domains &k,n and %, m are 
connected one with other by a chain of the pairwise intersecting domains ^k s ,n s - 

ii)% = {(k,0),keN p - 1 }. 

Hi) Ifr° k n < rl, n , and (k, n) G (k! , n') G %< , then j ^ f . 
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Then the clusters have the form 

<if j = {(2k m ,n m ),meN tj }, V-j = {{2k m - 1, n m ), m G for all jgN (4.12) 

where the number of elements of the cluster c ^'±j satisfies the estimate £±j ^ |, n m G N, 
1 ^ k\ < ... < kt ± , ^ |, 1 ^ rie ±j < ... < n\. Introduce the domains 



■ f .W kin , j g z. 



(4.13) 



These domains satisfy the following relations: 

% C {A G C: \lmz\ <P}, y~j C {A G C : | lmze~ i ^p \ < /?}, (4.14) 

for all j ^ 0. Moreover, if Xj G Aj/ G 'Vy and j < j', then ReAj < Re Xf, and each 
domain "Vj-\ is separated from "Yj,j G Z, by the line {A : Rez = Rj} for some Rj G R, 
... < R-i < R < Ri < R 2 < .... 

Lemma 4.3. The function p has as many zeros, counted with multiplicity, as the function p°, 
in each domain {A G C : R-n < Rez < Rn, \ lmz\ < N} and in each domain ~fj, \j\ > N for 
N G N large enough. There are no other zeros. 



Proof. Recall that 



P 



] | (cosh 



ZUj 



cosh zui)' 



2(p-i)p 

Assume that for each 1 ^ j < i ^ p, and for some c > 
I cosh zlo, - cosh zu t I > ce max{ 1 Re z ^ 1 ' 1 Re ZU)e 1 } as 



-£(£±1) „ x „_l OO p-1 



n=l fe=l 



fc,n 



^1 ^ fl 

z — 7m — > p, 



Asymptotics (I1.12p and estimates (I4.15P yield 



all (jfe, n) G N p _i x N. 

(4.15) 



P(A) 
P°(A) 



n 

1<7<^<P 



A,(A)-A,(A) 



cosh zu;^ — cosh zu>£ / 

coshz^ - cosh zw/ + 0(|z| _1 el Reza ^) + 0(|z| _1 el Re ^l 
cosh zujj — cosh 



0(\z\ 



as I A I 



n 

(4.16) 

s- 00. Let N G N be large enough and let N' > N be another integer. Let A belong to 

d'fn, |n| > A. Asymptotics (14 . 1 6 j) on all contours yields 

P(A) 



the contours Co(i?jv), Co(Rn> 

\p(X)-p (X)\ = p°(X) 



P°(A) 



1 



|p°(A)|0(|^r 1 ) < |AA)|- 



Hence, by Rouche's theorem, p has as many zeros, as p° in each of the bounded domains 
and the remaining unbounded domain. Since Ni > N can be chosen arbitrarily large, the 
statement follows. 

> p for all (k, n) G N p _i x N. Then 
Up+k+i) i i2im\ > 2(3ck and a fortiori 



2ii 



We have to prove estimates (I4.15p . Let \z — nn 
\2zrj k Ck — 2ixn\ > 2ftcp.. Identities (12.31) give \z(u p+ k 
\z(u>j ± Ue) + i2ixn\ > 2/3c/c for all 1 ^ j < I ^ p. Using the standard estimates we obtain 



cosh zujj — cosh zui \ = 2 



sinh 



z(Uj - Ut) 



sinh 



z{ojj + ue) 
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for some c > 0, which yields (I4.15p . ■ 

Lemma 4.4. i) Let r p+k (X) = r p+ j(X) for some 0^fc<j^p,AG A^, where R > is large 
enough. Then j — k + 1 and A G ^k,n for some (large) nGN. Moreover, in this case A G ^4 jn 
and A is also the zero of the function r p+k — r p+k+1 . 

ii) Let r£' , k G N p _i, n ^ 0, be ramifications of the Lyapunov function for the operator . 
Then 

<n =(-l) fc (^) 2P (l + (-l) P+1 ^(l+0(^ 2 ))) as n^oo, all k G N^. (4.17) 
Moreover, all ramifications rt 1 ' are real for n large enough. 



Proof, i) Let r p+k (X) = r p+ j(X) for some 0^k<j^p, Xe Ar. Then, due to Lemma fl~Tj 
j = k + 1. Asymptotics (11.81) gives 

e z(Q P+k -Q p+k+1 ) = i + o(\ z \~ 1 ) as |A|->-oo, where z = A^. (4.18) 

Substituting identity (12. 3 h into (14. 181) we obtain z = nn— + 0(n~ 1 ) as n — > oo. Then A G ^4 n 
for n G N large enough. 

The domain & k ,n is symmetric with respect to the real axis, then A G &k,n- The function 
p is real on R (see Theorem 11.21 i). then A is the ramification. Note that A is a zero of the 
function r p+k — r p+k+1 . Then asymptotics (14.11) and the first identity in (I2.13P show that A is 
also a zero of the function r p+k — r p+k+1 . 

ii) Let A = r^f n for some (k,n) G N p _i x N. By Lemma SHi), z = {r k fj^ = z + e, 
where z = (rl 2n )^p and \e\ < p(3 for all n G N large enough. Then 1 = r^ 1 (A)(r^ +1 ) _1 (A) = 
e «(n^(A)-n^ +1 (A))^ w h ere s = p _)_ ^ an( j we usec i identities (I2.2ip . Consider the case ImA ^ 0, 
the proof for the other case is similar. Then z satisfies the identity 

2nni 

Z = ^(A)-^ +1 (A)' (419) 

Asymptotics (I2.22j) yields 

z= + 0{n- 2 ) = ± '- + 0{n- 2 ) as n -+ oo, (4.20) 

w s — (jJ s +i c k 

where we used (I2.3p . Asymptotics (I2.22p and (I4.20p give 



<(A) - ^ +1 (A) = u . - - - J- | + 0(n 



T) 



, +1 W 1+ (l)^\ + 4 



where we used (12.31) and the simple identity w s w s+ i = r] k 2 . Substituting asymptotics (14.211) 
into identity (I4.19P we obtain 



d _ t*M\ + (Din-*) 



which yields (I4.17p . 
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Assume that r%' n are non-real for some n G N large enough. Then r%' n = r£'*, which is in 
contradiction with asymptotics (I4.17p . Hence e M. ■ 
We introduce the labeling of the ramifications at high energy: 

For each k G N p _i,n ^ uq for some uq G N large enough, r^ n are zeros of the function 
Tp+k — Tp+fc+i and r^ n G Moreover, we assume that lmr^ n ^ and 

l f lmr k,n > > theU r k,n = r k,W 

if^tn = 0, then (-l) fc r fen < (-l) fe r+ n . 
Corollary 4.5. T/ie following identities hold true 

T P+k(r% n ) = T p+k+ i(r± n ), r p _ fc (rjj = r p _ fc+1 (rjj, A p _ fe (r£j = A p _ fc+1 (rjj, (4.22) 
/or a// k G N p _i, n ^ no for some no G N /arge enough. 

Proof. The results of Lemmas 14.31 and 14.41 i) yield the first identities in (14.221) . Identities 
(14. 2 p give the second identities in (I4.22p . The definition of the functions Aj, see Theorem ll.il 
i), implies the third identities in (I4.22p . ■ 

Remark. Identities (I4.22p define the order of attachment of the sheets of the Riemann surface 
& at high energy (see Fig. [1]). 

5. Asymptotics 

Now we will determine the rough asymptotics of the periodic and antiperiodic eigenvalues 
and ramifications of the Lyapunov function. 

Lemma 5.1. The periodic and antiperiodic eigenvalues A^ and the ramifications r^ n satisfy: 

A± = (nn) 2p + 0(n 2p - 2 ) as n ->■ oo, (5.1) 

rt = rl n + 0(n 2p - 2 ) as n -+ oo, k G N p _i, r° n = (-1)*^* (5.2) 

Proof. Let A = A,^ for some n G N. Lemma 14.21 gives z = (A^)^p = nn + 5, where \5\ < | 
for all n ^ 1 large enough. The periodic and antiperiodic eigenvalues are real zeros of the 
functions — l,j G N p . Asymptotics ( 14. ip show that these functions, with only exception 
Tp — 1, have no any large real zeros. Then 1 = t 2 (A) = e 2lz (l + 0(n -1 )), where we used (14. ip . 
Substituting z = 7m + 8 into this identity we obtain e 2lS = 1 + Oin^ 1 ). Then 5 = 0(n _1 ) and 
z = irn + 0(n _1 ), which yields (15. ip . 
We will prove (15.21) . Let A = r^ 2n for some (k,n) G N p _i X N. By Lemma \A.2\ z = 

( r fc2n)^ = z o + £ ) where zq = (r^ 2n .)^ an d \s\ < p/3 for all n G N large enough. Identities 
(I4^2|) show that 1 = r s (A)r s+ \(A) = e^-^+i) (1 + 0(n^)), where s = p + k and we used 
asymptotics (14. ip . Substituting z = zq + s into this identity and using e z o(n s -n s+ i) _ ^ we 
obtain e <to a -n s +i) = \ + Oin" 1 ). Then e = O^" 1 ) and z = z + O^" 1 ), which yields ( l5T2"j) . 
■ 

In order to improve asymptotics (15.11) . (15. 2p we determine the asymptotics of the function 
D(t, A) = det(A4(l, A) — rl 2p ) in the neighborhoods of the unperturbed ramifications at high 
energy. 
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Lemma 5.2. Let k E and let X E C + and r E C satisfy 

A = r° n + CV^ 2 ), r = T p ° +fc (A)(l + 0(n- 1 )) as (5.3) 

Then 

T = r° +k+1 (\)(l + 0(n- 1 )) as n^oo. (5.4) 
Moreover, the determinant D(t, X) , given by ( [1. 7| ), satisfies the asymptotics 



D(t, A) = ct(r, A) det 



1 " W*(A))- x t 



fk,n 

1 " (**fi(A))- 



0(n 



(5.5) 



as n — >• oo ; where a = t p k 1 rij=p+fc+2(' r j 1 (A)) 1 7^ anc ^ fk,n are given by A3. 35) . 



Proof. Identities fl2T2TD give ftjf(A) = for all j G N 2p . We have z = X^ = z° + O^ 1 ), 
where z° = (r^ n )^ = vrn^-. Asymptotics (15. 3p and identities ( 12. 3 p give 

re -**+i = e*(«--".+0(i + o( n -i)) = e z °(^-^+0(i + o( n -i)) = l + O^ 1 ), s = p + A; (5.6) 

as n — )• oo. Asymptotics ( 15. 6 p yields (15 .4p . 
We will prove asymptotics (15. 5p . Identity ( I3.30P yields 



D(r, ■) = detOFe**" - rl 2p ) = det(.F - re - **") = det 
where the matrices Bi = diag(o;f, 5 2 = diag(a;^ +2 , ...,o;^,), 



A 3 .-. i 



A 2 

A 4 - re- zB2 



, (5.7) 



A x 



A, 







•F. s+i,i •• 


■ J~ S+l.S+l , 


r s+2,1 


•• J 7 s+2,s+] 




J 7 2p,s+2 



An 



A, 



-> 1,3+2 


J 7 l,2p 1 


J 7 s+l,a+2 •• 


•F. s+l,2pj 


'J 7 s+2,s+2 


■■ J 7 s+2,2p 


V J 7 2p,s+2 


J 7 2p,2p 



Dir. ■) — a det 



a = a(r, A). 



Due to (13.321) . the matrices Ai,...,Aa are bounded for |A| > large enough. Identity (15.71) 
yields 

f A x - re- zBl A 2 

Estimates (p?7j) . asymptotics (P2T25]) and (Q yield |r| -1 e* BaW = e - 2 ^e 2 ^+ 2 (l + O^ 1 )) as 
n -> oo. Relations <^M> show that Rez(u s -u s+2 ) > a\z\, a > 0. Then |r| _1 e* Ba W = 0(e~ an ). 
These asymptotics show that the matrix r~ l e 2 Aa — l p _fc_i is invertible for n large enough. 
Using the standard formula (see [Gaj . Ch.2.5) we obtain 

D(t, •) = a det (A 1 - re^ 1 - ^(^V* 2 A 4 - Vfc-i)" 1 ^"^* 58 ^) det (t~V B2 A4 - V*-i) • 

(5.8) 

Substituting the asymptotics |r| 1 e zB2 ( A ) — 0(e an ) into identity ( 15. 8 p we obtain 

D(T,X) = a(det{A l {X)-Te- zBlW ) + 0{e- an )) as n -+ oo, a > 0. (5.9) 
Furthermore, we have 

'As - re- zB * A 6 \ 

A 7 An - Te~ zBo J ' 



A x - re 



(5.10) 
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where B 3 = diag(e xu i,...,e 2a, «- 1 ), B = diag(e~ 

A* 





•F S,l ••• 3~ S,S— 1 1 A I ^~ 9,3 J~ S.S+l 




s-1, 




/I7 IX" 77 I ) I 77 77 

\' / ~ S + 1,1 ••• S + 1,S-1/ \*' S+1,8 S+l.S+l 

Estimates ([277]) imply re - zB 3(A) = ^-^.-i^ + = e *w.+i e -*w,_i( 1 + w h e re we 

used dS3D- Relations ((2HD show that Re^(w s _i - w s+1 ) > o|z|. Then T e~ zB ^ x) = 0(e~ an ). 
Asymptotics (13.321) show that ^}j(A) = 1 + 0(n~ 2 ), j G N 2p , which yields 

A 5 (X) - T e- zRi(x) = l s _i + 0(n" 2 ) as n -> 00. (5.11) 

Thus the matrix A 5 (X) — re^ zBz ^ is invertible for large n and f)5. lOj) gives 

det(Ai - Te~ zBl ) = det{A 5 - Te~ zB3 ) det(A - Te~ zB ° - A 7 (A 5 - Te~ zB:i )- 1 A^). (5.12) 

Substituting asymptotics (15.111) into identity (15.121) we get 

detiA^X) - re- zBlW ) = det(A (X) - re- zBoW + 0(n- 2 )) as n ->• 00. 

Substituting this asymptotics into (15.91) we have 

n( \\ a A ' (?s,s(X)-Te- z ^ ^,s + i(A) \,nf-^\ 

Substituting (I3.32p . (13.341) into the last asymptotics we obtain (15. 5p . ■ 

Below we write a n = b n + £ 2 (n) iff the sequence (a n — b n ) n ^i G £ 2 . 
Proof of Theorem 11.21 i) Repeating the arguments from [CKJ we obtain that p is entire 
and real on R. Asymptotics ( 11.12)) yields (I1.17p . 

ii) Let A = r^ n for some (k, n) G N p _i x N. We assume that Im A ^ 0. Then Qj = Uj and 

= for all j G N 2p . Using the identity r A T n = r^ n we obtain asymptotics for r~ kn G C_. 
Let A M = r£n, where \x = q Pj o, be the unperturbed ramification. Asymptotics (14.17j) . (15. 2p 
yield 

z = = z M + £<f, where z M = (A M )^, £ = — , 5 = OfrT 1 ) as n 00. (5.13) 

Cfc 

Since A = r^ n is the ramification, the monodromy matrix A4(1,A) has the eigenvalue r of 
multiplicity 2, i.e. its characteristic polynomial D(-, A) has the zero r of multiplicity 2. If n G N 
is large enough, then identities (I4.22p show that r = r p+ fc(A) = r p+ k+i(X). Using asymptotics 
(14. ip we obtain r = T° +fc (A)(l + Oin^ 1 )) as n — > 00. Then we can apply asymptotics (I5.5p . 

Note that the function Air) = det ( fll ™ 2 ' ' ) , where ou G C for all 7 G Nk, has the 

\ 0.4 a 5 — ra 6 I ' J jo, 

zero of multiplicity 2 iff (a 2 a 5 — ai<2 6 ) 2 + 4a 2 a 6 a 3 a4 = 0. Using asymptotics (15. 5p we deduce 
that 



(tfCA))- 1 ^ + 0(n- 2 )) - (rf +1 (A)) _1 (l + 0(n- 2 ))^ 

+ 4(^(A)< +1 (A))- 1 (/ M + O(n- 2 ))(7 M + O(n- 2 ))=0 as n -> 00, (5.14) 
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where s = p + k and rj 1 are given by (I2.2ip . Identity (I5.14p yields 

W(A))-*WVi(A)) 4 (1 + 0(n- 2 )) - (t»(\))H< +1 W)^ (1 + 0(n- 2 ))) 2 + 4|/ M | 2 = n- 3 £ 2 (n) 

(5.15) 

as n — > oo, where we used fk, n = i 2 (n)0(n x ), see (13.351) . 

Identities (14.221) . applied to the operator H^, yields r s M +1 (A M ) = r s M (A M ). Asymptotics (15.131) 
imply 

^(A)« +1 (A))- 1 = ^{\){^ +1 {\))-\r^\n)- l r^) 

= e ^(^(A)^(AM)-^ +1 (A)+< +1 (AM)) e ^(A)-^ +1 (A))^ /g 16 x 

Asymptotics fl2~25|) gives ojf( A) = w j (A) + 0(n- 2 ). Asymptotics fl2~23|) yields w?( A) -otfiX 1 *) = 
0(n~ 4 ) as n — > oo, j G N2 P . Then (15.161) implies 

Ctf (A))* = effKW- w -+iW)(l + 0(ri~ 2 )) = 1 + «(-l) fc+1 <5 + 0(n~ 2 ) 

as n — > oo, where we used (I2.3p . Substituting these asymptotics into (I5.15P we obtain 

5 = ±(-l) k \f kn \ + 0(n- 2 ) = ± ^ l ) k °^\ + Q ( n -2) as n ^ 00; 

2p7T7i 

where we used (I3.35p . Then 

r± ^ = «a* + t- 

2p7m 



and 



± 1 *'" V X; ' /fcl ^"' +0(^- 3 )=r^±( — ) c fc |g p , n |+0(^- 3 ), (5.17) 

where we used z^ = £7rn + 0(n^ 1 ), see (I4.17p . and % p = (— l) k , see (12. 4p . Substituting 
asymptotics fl4~TTj) into (Q7jl we obtain Q| . ■ 

Proof of Theorem 11.31 i) Asymptotics f II . 1 2 j) shows that the branches of the Lyapunov 
function A on the interval (K, +oo) for some K G M satisfy: 

if p is odd, then there is exactly one real branch Ai and the other branches are non-real; 

if p is even, then there are two real branches Ai and Ae +1) Ae +1 > 1 and the other branches 
are non-real. 

Moreover, A X (A) ^ Aj(A) for all j = 2,...,p, A G (K,+oo), hence Ai is analytic on the 
interval (K, +oo). Asymptotics (I1.12p for A x and Theorem 11.11 ii) show that the function Ai 
oscillates on (K, +oo) similar to the Lyapunov function for the Hill operator. Then using 
identity (I1.13P and the standard arguments (see |BKlj ) we obtain the needed statement. 

ii) Let A = A^ for some n G N. Recall that A G 1R and satisfies 

D n = D((-l) n , A) = det(A4(l, A) - (-l) n l 2p ) = 0. (5.18) 

Let A M = \% ,ziz be the unperturbed 2-periodic eigenvalues. Asymptotics (12.241) . (15.11) give 
z = = z^ + 5, where z^ = (A M )^, 5 G R, 5 = C^n" 1 ) as n ^ oo, A M satisfies (I2.24p . 



M)2p-i(_i)^ fc | $ | /7m\ 2p - 2 
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Asymptotics (15. 5p for k = give 

A. = a.d rt (( 1 -(- 1 £W A »" 1 _,_ 1) ^, +l(A)) -.)+0(n-)) a s 

(5.19) 

where a n = «((— l) n , Aj) 7^ 0, /o, n = — ^J^. Each A^ 1 is a periodic or antiperiodic eigenvalue 
of the operator and the corresponding multipliers satisfy the identities t^A^ 1 ) = Tp +1 (A M ) = 
(-l) n . Then 

(-l) n (rj i (A)) _1 = (rj t (A))- 1 r;(A^ t ) = e -(^«tfW e *MC*0 = e -K e -^KW-^(A")) ( 5 . 20 ) 

for j = p,p+l. Asymptotics (I2.25j) yield e~ Sw ^ x> = e~ SuJj (1 + 0(n~ 3 )) as n — >■ 00. Asymptotics 
(12^3]) show that wj*(A) - w?(A") = 0(n" 4 ). Then flCTj) gives 

( T M(A))- 1 = (-l)" e -^(l + 0(n- 3 )) = (-l) n (l-5a; i )+0(n- 2 ) as n ->• 00, j = p,p + l. 
Using the identities u p = — u p+ i = —i we obtain 

(T p u (\))~ 1 = (-l) n (l + iS)+0(n- 2 ), (r^ 1 (X))- 1 = (-l) n (l-i5) + 0(n- 2 ) as n -)• 00. 
Substituting these asymptotics into (I5.19P we obtain 

D n = a n det((^ 6 fo A + 0(n~ 2 )] = a n det(Af + 61 2 + 0( n - 2 )), Af = ( ^ 

\\Jo,n ld J J V l J0,n 

Using identity (I5.18P we conclude that 5 is an eigenvalue of the matrix M + 0(n~ 2 ). Since 
the eigenvalues of the matrix M have the form ±|/ 0)n |, we obtain 5 = ±|/o,n| + 0(n~ 2 ). Using 
the identity \f 0>n \ = ^ we obtain \& = z>* ± + 0(n- 2 ). Then (|2"23|) gives (PUD . ■ 
Proof of Corollary 11.41 i) Asymptotics (ll.20p and the estimates \g p>nh \ ^ give the 
asymptotics 

| 7n | = A+ - A; = 2(7m) 2p - 2 \q p , n \ + 0(n 2p - 3 ) as n -> 00, 

and the estimates |7„J ^ 27r 2p_2 n^ p_2_a (H-0(n/ c a ~ 1 )) as fc — > 00, which yields the statement. 

ii) Assume that r^ n . are non-real. Then r A T n . = r^ n . , which is in contradiction with asymp- 
totics f ll.lSp . Hence rf n . G R. Moreover, asymptotics fll.ISp gives 

c k\q P ,n\ + 0(n 2p_3 ) as n — > 00, 

which yields — \ — > 00 as k — > 00. ■ 

6. Appendix 

Proof of Lemma 13.21 The standard arguments yield that the fundamental solution A4(t, A) 
of equation (11.50 . with the initial condition A4(0, A) = l 2p , satisfies the integral equation 

M(t, A) = Mo(t, A) + / .M (*-s,A)Q(s).M(s,A)ds, (6.1) 

where -Mo(^, A) = e* P( ^ is a solution at Q = 0. 

Describe the properties of the matrix-valued function .Mo- Each function A4o{t, ■), t G M, is 
entire. Moreover, A4 = (( < / ? j)^ fc_1 ^)fe^=i> where </?°,j G N 2p , are the solutions of the equation 



k,rij 



2(1 

Cu J 
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y( 2 P) = satisfying the conditions (cpj) ( k ^(0, A) = 5jk for all k G N 2p . These solutions are 
given by the identities 

1 JO 

Then K^)^- 1 )^, A) | < e 20 l*l for all j, fc G N 2p , (t, A) G E x C, and 

\M (t, A) | < 2pe zo141 , all (t, A) G R x C. (6.2) 

Estimate (16.21) will be useful for obtaining the first estimate in (13. lip . Now we will prove the 
other estimate of A4q (see below (16.31) ). which will be effective to obtain the second estimate in 
(13. lip . In fact, direct calculations show that V = zZCB{ZC)~ 1 , where the diagonal matrix B 
is given by B = diag(uj)l P , and the matrix C has the form C = (u)*~ l )^ =v Then C^ 1 = j^C*, 
Mo = ZCe ztB (ZC)- x and 

\Z-\X)M (t, \)Z{\)\ 2pe z °W, all (t, A) G R x (C \ {0}). (6.3) 

The standard iterations in (16. ip yield 

M(t,X) = Y.M n (t,X), M n (t,X)= f M (t-s,X)Q(s)M n - l (s,X)ds. (6.4) 
and A4 n (t, A) is given by 



M n (t,X) 



J n^Q(Wi-*fc,A)Q(t fe )^ Q (t 1 ,A)dt 1 dt 2 ...^ n (6.5) 



the factors are ordering from right to left, T = {0 < t\ < ... < t n < t n+ i = £}. Substituting 
estimates (16.21) into identities (16. 5p we obtain 

\M n (t,X)\ ^ -^ eZ ° ltl (iP f Q \Q( s )\ ds ^j , a11 (ji,t,X) G N x R x C. (6.6) 

These estimates show that for each fixed i 6 1 the formal series (16. 4p converges absolutely 
and uniformly on bounded subset of C. Each term of this series is an entire function. Hence 
the sum is an entire function. Summing the majorants we get 

\M(t,X)\ ^ 2pe zoltl+ ti lQ{s)lds all (t,X) G R x C, (6.7) 

which yields the first estimate in (13. lip . Moreover, we deduce that the monodromy matrix 
A) is a continuous function of all qj G L 1 (T), j G N p . 
Substituting estimates ( 16. 3p into identities (I6.5P we obtain 

\Z- 1 (X)M n (t,X)Z(X)\ *C ^e^Up j\z-\X)Q{s)Z{X)\ds 

for all (n, t, X) G NxRx(C\{0}). Substituting these estimates into the series (16 .4p we obtain 

\Z- 1 (X)M(t,X)Z(X)\ ^ 2 pe ZoW+ ti lz ~ 1WQ{s)zwlds , all (t, A) G R x (C\{0», 

which yields the second estimate in (13.111) . ■ 
Proof of Lemma 13.41 Asymptotics ( I2.25P implies 

U(X) = Z(X)U (l 2p + 0(\z\- 2 )), where Z = diag(^- 1 )f , U = (W k 7 l )% =1 . 
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Then 

Q(£, A) = W _1 (A)Q M (£, X)U(X) = U a l Z- l {\)Q^{t, X)Z(X)U (l 2p + 0(|z|- 2 )) as |A| oo 
uniformly on i G [0, 1]. We have 

Z-\\)Q?(t,\)Z(\) = ^ — [(q p (t) -n)E + b(t)0(\z\- l )j as |A| oo, 

uniformly on t G [0, 1], where S is given by ( ESP - Then we obtain f l3T8|) with £ = W _1 £W 
Using the identity Wq" 1 = we get C = ^UqEUq. Then 

2 P nfnP-i 



which yields the identity for £ jfe in (|3.18p . ■ 

Proof of Lemma 13.51 i) Asymptotics (I3.18p show that the estimate 

max \Q kj (t,X)\ ^ maX ^^ ( ^ p^llll , all M)ex[0,l]xA Rl (6.8) 

for some i?i > 0. Assume that 

max |ejy(t,A)| ^ 2, all (£, A) G K x A^, for some i? 2 > 0. (6.9) 

fcjeN 2p 

Then substituting estimates (16. 9p into (I3.23P we obtain 

2p r 1 2p 

\\(KQA)(-, X)\\oo = max \(KQA){t, X)\ = max V/ e^t-s, X) V] Q iri (s, X)A nj (s)ds 

2p r 1 ~ 2p r 1 ~ 

^ 2 max V / |Q in (s,A)||Ay(s)|ds < 2 max V max \A nj (t)\ / |Q in (s, A)|ds. (6.10) 

i,n=l i,n=l 

for all A G A max { fil) ^ 2 }. Using estimates (16.81) we obtain 

c 2P 
\\(KQA)(-, A)||oo < ■pr max max V|Ay(£)|, 



n=l 



which yields (13. 26 p . 

We will prove (16. 9p . We will consider the case % > j. The proof for i ^ j is similar. Identities 
(I3.24p shows that ey(£, A) = for £ < 0. Asymptotics (I2.25P shows that 

= fL + — , 7 G N 2 «, where max l^-fA)! < - log 2 

for some R 2 > 0. Estimates (12.71) give Rez(f2j — fl,) ^ for i > j. Then 

(^)| = e iRe Z (nf(A)-^(A)) ^ e me z{Oi-Oj)+t Re^A)-^ (A)) ^ e |^(A)-^(A)| ^ 2) (6.11) 



for £ ^ 0, i > j, \z\ > R 2 , which yields estimates (16. 9p for % > j. 
ii) The standard iterations in (I3.22p give 

oo 

G = Y<Gn, £o = l 2p , G n = KQg n _ 1 = (KQ) n , all n G N. (6.12) 

n=0 
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Estimate ( jSjZgJ) yields 

||^(-,A)|| 00 = ||(^Q)"(-,A)|| 00 ^ (JL^\ , all (n,X)ENxA max{RM . (6.13) 

These estimates show that the formal series (16.121) converges absolutely and uniformly on any 
bounded subset of A R ,R = max{i? 1; R 2 , £ 2p }. Hence it gives the unique solution of equation 
(I3.22p . Each term of this series is analytic in A^. Hence the matrix- valued function Q is 
analytic in A^. If A G A 2 _r, then substituting estimates (16.131) into the series (16.121) we obtain 



n=0 I / 1 \z\ n=l ^ I I ' 

°° f P \ n { — ) 2 9P 

\\G(; A) - l 2p - G 1 (; A)|U ^ (±r) = ^L- ^ ,A 



n=2 v 1 1 7 |*| 

which yields (I3.27p . Substituting (13. 18j) into the identity Q\ = KQ we obtain 



(^i(t,A))y = ^ / e ij (t-s,X)(q p (s) - fi)C tj (X)ds + a^t, X)0(\z\- 2 ), i,j G N 



2p> 



where aij(t,X) = (— l) p+1 J ey(t — s, X)b(s)ds,i, j G N 2p . Estimate (16. lip gives |Oy(t, A)| ^ 
2||6||, which jdelds 

iii) The matrix-valued function Q(t,-),t G [0, 1], is analytic in A^, then T(t, •) is also analytic. 
The second estimate in (I3.27P yields asymptotics (I3.29p . ■ 
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